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<L> Abstract 



In this work, using self-consistent tight-binding calculations, for the first time, we show that a direct 
to indirect bandgap transition is possible in an semiconducting armchair graphene nanoribbon by the 
application of an external bias along the width of the ribbon, opening up the possibility of new device 
applications. With the help of Dirac equation, we qualitatively explain this bandgap transition using the 
asymmetry in the spatial distribution of the perturbation potential produced inside the nanoribbon by 

i— I 

J> the external bias. This is followed by the verification of the bandgap trends with a numerical technique 

ON 

t-H using Magnus expansion of matrix exponentials. Finally, we show that the carrier effective masses possess 

tunable sharp characters in the vicinity of the bandgap transition points. 
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Graphene, the two dimensional allotrope of carbon, has drawn an enormous amount of at- 
tention in the literature after its first isolation on an oxide substrate [JJ [2]. Apart from the 
theoretical physics point of view [3l Uj , graphene has emerged as a possible candidate for dif- 
ferent electronic devices including field effect transistors [2]-[B]. However, the small bandgap of 
graphene reduces the controllability of such devices and thus limits its widespread applications. 
Graphene Nanoribbon (GNR), on the other hand, a quasi-one dimensional strip of graphene, 
has been shown to provide a significant bandgap [9]- [12] and hence is being considered as a 
channel material in field effect transistors |13j-|15|. 

More recently, it has been theoretically shown that the bandgap of a graphene nanoribbon 
can be tuned significantly by the application of an external field along the width of the ribbon 
[TBI ITT1 [TH] and at sufficiently large field, it is also possible to collapse the gap of the ribbon. 
In this work, we extend this result in a more generic way for both semiconducting and metallic 
Armchair Graphene NanoRibbon (A-GNR) and demonstrate that not only the magnitude of 
the bandgap, but the whole electronic structure of the nanoribbon can be altered significantly 
depending on the magnitude and polarity of the external bias. In particular, we will show 
that it is possible to obtain a bias dependent direct to indirect bandgap transition in such a 
nanoribbon. This kind of control over the electronic structure by the application of an external 
bias opens up the possibility of new device applications. 

A schematic of the setup that we consider in this work is shown in Fig. [TJ where an Armchair 
Graphene NanoRibbon A-GNR of width W is sandwiched between a left gate (G^) and a right 
gate {Git). There is a third contact V c which keeps the chemical potential of the A-GNR at 
zero. In this work, we will primarily focus on semiconducting nanoribbons, i.e., the number of 
dimers (N) along the width, is of the form 3M and 3M + 1. Note that, N = 3M — 1 gives 
rise to metallic nanoribbons [T2"] and is briefly analyzed in this work. The gate terminal in 
each gate stack is separated from the GNR by a dielectric. We assume the Equivalent Oxide 
Thickness (EOT) of the gate dielectric to be lnm. The interfaces between the GNR and the 
gate dielectric are assumed to be perfect, hence the electronic structure of the GNR is not 
altered significantly. The gate dielectric confines the electrons and holes in the GNR along the 
x direction, however the y component of the states can be obtained using longitudinal wave 
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vector k y = k. The work function of the gate metal is assumed to be such that a zero hatband 
voltage is obtained. The external biases Vi and V r at the left and the right gates respectively 
can be varied independently. The application of an external bias changes the potential energy 
U(x) = —q<f>{x) inside the GNR, altering the electronic structure. 

We now provide the details of the calculation procedure used to investigate the effect of 
external bias on an A-GNR. The self-consistent electronic structure of the A-GNR is determined 
by using tight binding method ([19, 20J) coupled with the Poisson equation. Taking the left 
edge of the GNR at x = and the plane of the GNR as z = 0, the charge density is given 
by p(x,z) = qn(x,z), where q is the electronic charge and n(x,z) is obtained as the difference 
between the hole [rih(x,z)] and electron [n e (x, z)\ density as 



n(x, z) 



^(l-/(^(fc)))|^(x,z)| 2 -^/(^(fe))|^|(x,z)p 



I 2 (1) 

i,k j,k 

where f(E) = l+e (E-^)/k B r ls the Fermi-Dirac probability at temperature T. Here k goes over 
the whole first Brillouin Zone, i and j are the valence and conduction band indices respectively. 
The chemical potential p, set by the contact V c , is taken to zero. Eiik) is the energy eigenvalue 
of the state (i,k) obtained from the tight binding bandstructure taking only p z orbital into 
account, with an intra-layer overlap integral, S = 0.129 between two nearest carbon atoms 
and the intra-layer hopping t as — 3.033eV pj)|. Note that, the results obtained from the 
nearest neighbor calculation are in close agreement with simulation that take into account 
coupling terms up to the third nearest neighbor (see supporting information). To obtain the 
wavefunction ip£(x, z), we assume normalized Gaussian orbital as the basis function, where the 
parameter of the basis function is fitted using the parameter S. The wavefunctions are set to 
zero at the dielectric interfaces indicating an infinite potential barrier. Once self-consistency 
is achieved between the bandstructure calculation and the Poisson equation for a given gate 
bias, the energy eigenvalues at different k points correspond to the electronic structure of the 
A-GNR. 

We take an A-GNR with = 36 (W = 4.55nm) and consider three representative bias 
conditions, namely, (i) Vi = —V r , (ii) V\ > 0, V r = and (iii) Vi < 0, V r = 0. We now present 
the results in these three cases as shown in Fig. [2]|4| 



Case (i): In this case, the two gate voltages are anti-symmetric in nature, i.e., V\ = —V r = 
V g . We observe a significant reduction of bandgap in Fig. [2^a) and (b) with an increase in V g , 
and this has also been predicted in [T71 [18]. It is observed that at sufficiently large V g , both 
the conduction and valence band edges shift from k = 0, giving a 'Mexican Hat' shape around 
k = 0. Fig. [2^b) clearly shows a threshold like behavior of the bandgap change |17j . and as 
the band edges shift from k = (non-zero A/c), the bandgap starts decreasing significantly 
with bias. However, the particle-hole symmetry is almost conserved (the small asymmetry 
in Fig. [2^a) is due to the non-zero overlap S assumed between two nearest neighbor carbon 
atoms in the honeycomb lattice) and hence the bandgap continues to remain direct in nature, 
at any bias condition. Note that, the anti-symmetric bias condition forces the GNR to retain 
its charge neutral condition with similar electron and hole density, keeping the total effective 
charge density very low. 4>(x), dictated by the Poisson's equation, thus remains almost linear 
(uniform field) along x, as shown in Fig. [2^c). Note that bias dependent bandgaps match very 
well with one of the previously published reports based on non-selfconsistent calculations |18j 
and this linearity of (j)(x) along the width of the nanoribbon is the reason of this unexpected 
close match. 

Case (ii): In this case, the left gate is kept at positive bias, keeping the right gate grounded 
and the results are shown in Fig. [3^a)-(c). The bandstructure shows a dramatic change as 
compared with case (i). With an increase in V/ = V g , the conduction band minimum shifts 
from k = 0, giving rise to a 'Mexican Hat' shape around k = 0. However, the valence band 
maximum continues to remain at k = 0, irrespective of V g . Thus, at any V g for which Ak > 0, 
the GNR has an indirect bandgap. The direct and indirect bandgap regions are indicated 
in Fig. [3](b). The magnitude of the bandgap continues to show a threshold-like behavior as 
before, with an increased sensitivity of bandgap when the system becomes indirect. Note that, 
the spatial distribution of 4>{x) along the width of the nanoribbon is severely non-linear (non- 
uniform field) to support the increased electron density inside the GNR that arises due to the 
conduction band edge moving closer to the chemical potential. We will later point out that it 
is this strong non-linearity of 4>(x) that causes such a direct to indirect bandgap transition. 
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Case (iii): A similar case like (ii) can be constructed where the left gate is at negative bias, 
with the right gate grounded. The bandstructure in such a scenario is shown in Fig. [4jja) where 
the conduction band minimum continues to remain at k = 0, and the valence band maximum 
shifts away from k = 0, depending on the external bias. In this case, the valence band edge 
moves closer to the chemical potential resulting in a relative increase in the hole density. 

Note that the corrections due to the second and the third nearest neighbor interactions 
contribute at relatively large values of k, away from the zone center [21]. However, the band 
edge shift (Afe) from the zone center (k = 0), in all the above scenarios, are small (~ 5%) 
compared to the size of the Brillouin zone. Hence, the calculations with nearest neighbor 
interactions are accurate enough to predict such direct to indirect bandgap transition. 

In Fig. |4]^b), we generalize this result and show the transition from direct to indirect bandgap 
in the (Vi,V r ) space. We compute the absolute difference of the k values of the conduction band 
minimum and the valence band maximum for any arbitrary combination of Vi and V r . This is 
plotted as a function of (V h V r ) in Fig. gb). A zero value (dark color) indicates direct bandgap, 
whereas a non-zero value (lighter color) represents indirect bandgap region. We clearly observe 
that in the (Vi,V r ) space, there are symmetric pockets of indirect bandgap regions, with the 
chosen cases (ii and iii) are the most favorable conditions to obtain such a bandgap transition. 

In the case of a metallic A-GNR, it is interesting to note that an asymmetric external 
electric field along the width opens a small bandgap at the zone center. Fig. [5^a) shows a 
direct bandgap of ~ 18meV for a metallic A-GNR with A = 35 under a bias of 2.8V at the left 
gate, while grounding the other. However, as shown in Fig. [5^b), at larger bias, this bandgap 
tends to become indirect accompanied with a reduction in its magnitude. We do not observe 
such an effect in the case of anti-symmetric bias condition, in agreement with |17j . 

The external bias dependent direct to indirect bandgap transition, coupled with the change 
in magnitude of the bandgap can have significant effects in phenomena including band-to-band 
tunneling, electron-phonon interaction and optical properties. Such an external bias dependent 
tailoring of the electronic structure can provide us with the possibility of a wide variety of 
fascinating electronic and optoelectronic device applications. 
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Now, to get more insights, we present a theoretical analysis of the phenomenon by starting 
from the Dirac equation [T2| [T7] . We write the low energy states *&(f) = e lk ° x ipj t .{f) + 
e~ ° x ip-(r) in terms of smoothly varying envelop ip = ip-}- t P+ an d V'- have components 
on the A and B sublattices in the honeycomb lattice with ko = — 47r/3ao and ao = 2.44nm |17j . 
By making the replacement k x — > —id x in the Dirac Hamiltonian [12] , we can write Hip = Etjj 
where the Hamiltonian {H) for the nanoribbon is given as 

with H± = ±ihva x d x — hvka y — q(p(x)I. Here, a are the Pauli matrices and v ~ 10 6 m/s. 
To keep the analysis simple, we assume the intra-layer coupling parameter S to be zero. The 
armchair boundary condition with ideal edges forces 

V + (0)+V-(0) = (3) 

and 

7P + (W) + e ikoW TP-(W) = (4) 

Now, we give a simple argument to show why we observe a direct to indirect bandgap transition 
in setup (ii) and (iii), whereas setup (i) provides direct bandgap independent of external bias. 
If we write the full Hamiltonian H by discretization of space along x, we find, 

Tr{H) = -Aq^^i) (5) 

j 

which is equal to zero in case (i) and this holds good for any k. This is due to the anti-symmetric 
nature of the external bias and hence of (j)(x) about the mid point of the nanoribbon. Now, using 
the fact that the sum of the eigenvalues equals the trace of H, this condition forces the sum of 
the energy eigenvalues at any k to be zero. Thus the conduction band and valence band remain 
symmetric about forcing the bandgap to be direct at any external gate bias. However, in 
cases (ii) and (iii), the asymmetric gate biases introduce consequent asymmetry in the spatial 
distribution of (j)(x) and hence force Tr(H) to become nonzero, allowing asymmetry in the 
conduction and the valence bands. This manifests as a bandgap transition in the nanoribbon. 
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We now provide an independent numerical method derived from Eq. [2] to re-calculate the 
bias dependent electronic structure and verify the trend of direct to indirect bandgap transition 
obtained from tight binding calculations. To do this, We rewrite Hifi = Etfi as |17j 

d^ ± = ±C^± (6) 

where 

C(x) = ka z - ia x {q(p{x) + E)/hv (7) 

Since, in general, C( x ) does not commute for two different x, we can write the solutions in terms 
of Magnus series |22|I23|: 

^ ± (W) = e e± V±(0) (8) 

where 0± = Yl'jLii^V^j- @j ^ s the 3 t erm i n the Magnus series with the first three terms are 
given as 



n 



C,{x 1 )dxi, 



o 



rW rxi 
6 2 = dxi / dX2[((xi),((x2)}, 

Jo Jo 

rW pxi (-X2 

9 3 = dn dx 2 dx 3 ([C(x 1 ),[C(x 2 ),ax 3 )]} + [ax 3 ),[C(x 2 ),C(xi}}) (9) 
Jo Jo Jo 

Using Eqs. |3j |4]and[8j we obtain 

(e e + -e ikoW e e -)4; + (0) = (10) 
To get non-trivial solutions for ^+(0), we obtain 



det 



e e+ _ e ik w e e 



(11) 



For a given k, the set of values of E satisfying Eq. 11 gives the required energy eigenvalues, 
which can be found numerically. We have verified that the results obtained using this method 
show a direct to indirect bandgap transition in setup (ii) and (in) whereas the GNR continues 
to remain a direct bandgap semiconductor in case (i). 

As a special case, Q{x) commutes for two different x for k = and hence 9j(k=0) becomes 
zero for j > 1. Hence, we can readily observe from Eq. [lljthat as long as (p(x)dx = 0, we do 



not have any change in the energy eigenvalues at k = for any arbitrary (p(x). This is why we 
should not expect any change in E(k=0) for any external bias as long as V\ = —V r . Note that, 
in reality, as shown in Fig. [2^a), we do see a small change in E(k=0) under gate bias, which 
arises from non-zero overlap parameter S. However, in case (ii) and (iii), nonzero J <ft(x)dx 
introduces a bias dependent upward or downward shift in the E(k = 0) value depending on the 
polarity of the terminal bias. 

As a final comment, we extract the effective mass values at the conduction band minimum 
and the valence band maximum for case (i) and (ii) using the E — k relationship obtained 
from self-consistent tight binding calculations. The results are shown in Fig. [6](a)-(b). In both 
the cases, we observe strong non-monotonic behavior of the effective mass values, both for the 
electrons and the holes. In case (i) [Fig. [6^a)], the effective mass of the electrons follows that 
of the holes for both small and large biases. However, at some intermediate gate bias, where 
the band edges start shifting from k = 0, we notice a significant difference in the electron and 
the hole effective mass values, though both of them strongly peak about that point. A similar 
behavior is observed in the electron effective mass in case (ii) [Fig. [6^b)], which sharply peaks 
around the direct to indirect band transition point indicating a "flattening" of the conduction 
band edge when it moves away from k = 0. However, we do not observe any such sharp peak 
in the hole effective mass around the bandgap transition point. This sharp notch like behavior 
of the effective masses (with almost an order of magnitude change) is a unique feature of the 
influence of the external field on the electronic structure of A-GNR where one can selectively 
"slow down" the carriers by choosing the appropriate external bias condition. 

To conclude, using a self-consistent tight binding calculation, we have demonstrated that 
it is possible to change the bandgap of an semiconducting A-GNR from direct to indirect by 
adjusting the external biases at the left and the right gate, opening up the possibility of new 
device applications. Such a direct to indirect bandgap transition has been explained, both 
qualitatively as well as quantitatively, by starting from Dirac equation, to support the findings 
obtained from tight-binding calculations. Finally, the external bias dependent carrier effective 
masses have been shown to have non-monotonic sharp behavior around the direct to indirect 
bandgap transition point. 
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Figure 1: A schematic of the A-GNR device with the external bias along the width of the GNR. 
The chemical potential of the GNR is set through the contact V c which is kept at zero bias (the 
reference potential). The gate biases V\ and V r create an external field along x inside the GNR 
which tunes the electronic structure of the GNR. 
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Figure 2: The bandstructure of the A-GNR device (JV = 36) using tight-binding method with 
the external biases at Vi=— V r =V g =lAV. The Chemical potential, /jl, is set to zero in all the 
figures, (a): E-k relationship (in blue) in the vicinity of k=0 shows a significant reduction of 
bandgap, but the gap remains direct. The red lines show the case when V g =0. k is in terms 
of ^ where a = 3«o with ao = 1.42A (b): Change in bandgap E g as a function of V g shows a 
clear threshold-like behavior. In the same plot, we show the shift (Ak) of the band edges from 
k = 0. (c): (p(x) plotted along x inside the A-GNR shows an almost linear spatial variation, 
indicating a uniform electric field along the width of the ribbon. 
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Figure 3: The tight binding results of the A-GNR device with external bias conditions, 
V;=V^=2.8V, V r =0\. (a): E-k relationship (in blue) in the vicinity of k=0 shows that the 
conduction band minimum is shifted from k=0, while the valence band maximum is at k=0, 
causing an indirect bandgap. The red lines show the case when V g =0. (b): Change in the 
bandgap E g as a function of V g again shows a threshold-like behavior. In the same plot, we 
show the shift (Ak) of the conduction band minimum from k = 0. (c): <j){x) plotted along x in 
the A-GNR shows significant nonlinearity in it spatial distribution, indicating a non-uniform 
electric field. The Chemical potential /j, is set to zero in all the figures (a)-(c). 
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Figure 4: (a): E-k relationship (in blue) in the vicinity of k=0 for Vi=V g =— 2.8V, V^=0V 
shows that the conduction band minimum remains at k=0, though the valence band maximum 
shifted from k=0, causing an indirect bandgap. The red lines show the electronic bands for the 
unbiased condition, (b): Plot of the absolute of the difference between the k values (in units of 
2n/a) of the conduction band minimum and the valence band maximum in the (Vi,V r ) space. A 
zero value of the same indicates direct bandgap (dark color) , whereas non-zero values represent 
indirect bandgap (lighter color). We observe symmetric pockets of indirect bandgap regions in 
the (Vi,V r ) space. 
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Figure 5: (a): E-k relationship in the vicinity of k=0 (for Vi=V g =2.SV, V^-=0V) shows the 
opening of a small direct bandgap (~18meV) at the zone center, for a metallic A-GNR with 
N = 35. (b): The bandgap becomes indirect with a slight reduction in its magnitude (E g ~ 
14meV) atVi = 6V. 
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Figure 6: (a): Effective mass at the conduction band minimum and the valence band maximum 
as a function of external bias in case (i) with Vi = —V r = V g . The effective masses are 
strongly non-monotonic with the peak occurs in the vicinity of the bias point where the band 
edges shift away from k = 0. (b) Effective mass at the conduction band minimum in case (ii) 
(Vl = V g > 0, V r = 0) shows a sharp peak around direct to indirect bandgap transition point. 
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Figure 7: Supporting Information: The tight binding bandstructure of the semiconducting 
A-GNR with N = 36 under external bias conditions, V^=V^=2.8V, V r =0Y considering (a): 
nearest neighbor interaction, and (b): three nearest neighbor interactions. The E-k relation- 
ships (in blue) shows that the conduction band minimum is shifted from k=0, while the valence 
band maximum remains at k=0. The red lines show the bandstructure without any external 
bias. The Chemical potential is set to zero. The nearest neighbor calculation is found to 
accurately predict the direct to indirect bandgap transition. 
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